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Abstract 

In this work, a study of the average transverse momentum (p T ) as a function of the mid-rapidity charged hadron multiplicity 
(N c h) and hadron mass (in) in p-Pb and Pb-Pb collisions at LHC energies is presented. For the events producing low N ch , the 
average p T is found to scale with the reduced hadron mass, i.e., mass divided by the number of quark constituents (m/n q ), this 
scaling also holds for inelastic pp collisions at RHIC and LHC energies. The scaling is broken in high multiplicity p-Pb and Pb-Pb 
collisions, where, for (dlV c h/d7/) < 60 the average pj is higher for baryons than that for mesons, though they increase linearly with 
m/n q . This behavior is qualitatively well reproduced by Pythia 8, but not by hydro calculations, where an universal scaling with 
the hadron mass (and not with m/n q ) is predicted for all the multiplicity event classes. Only the central (0-60%) Pb-Pb collisions 
behave as expected from hydro. 

Keywords: Multi-parton interactions, color reconnection, flow-like behavior, nucleon-nucleon reactions, nucleon-nucleus 
collisions. 


1. Introduction 

The main goal of the heavy ion program is the study of matter under extreme conditions of high energy density, 
where QCD phase transitions characterized by deconfinement and chiral symmetry restoration (QGP), are expected 
to occur 00. To extract the genuine properties of the hot matter, the program includes studies in control experi¬ 
ments, like nucleon-nucleon and nucleon-nucleus collisions, where final state effects were not supposed to be present. 
Unexpectedly, recent results at the LHC reveal the presence of QGP-like effects, namely, long range angular correla¬ 
tions 0] and flow-like patterns J4j in the systems created in high multiplicity pp and p-Pb collisions at yfs — 7 TeV 
and V s 'NN = 5.02 TeV, respectively. 

Such features of data are qualitatively described by hydro calculations 0(6), where the formation of a medium 
and its subsequent fast thermalization are implicitly asummed, but not fully justified. For example, it is known that 
the validity of a hydrodynamic description requires small values of Knudsen number (<s 1), however, for the systems 
created in proton-nucleus collisions at sJsnn = 5.02 TeV, the Knudsen number is already close to one in almost all the 
space-time points for temperatures T < 165 Me V 0. Albeit, based on ideas of string percolation 0(9) it has been 
argued that deconfinement can also occur in small systems, the origin of the phenomenon has not been established ED- 
With the same level of accuracy as hydro, other calculations, which do not require the formation of a mini QGP, also 
do a good job to describe the data. For instance, color glass condensate (CGC), the QCD formulation on the limit of 
high occupancy number at low x-Bjorken, predicts the long range angular correlations and their azimuthal anisotropies 
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Figure 1. (Color online). Average transverse momentum as a function of the scaled hadron mass (m/n q ) for pp collisions at Vs = 7 TeV simulated 
with Pythia 8.201. Two multiplicity event classes are shown, low multiplicity (left) and high multiplicity (right). Results are presented for 
simulations with (full red markers) and without (empty blue markers) color reconnection. The solid (dashed) line indicates a linear fit to meson 
(baryon) data. 


for even harmonic^] and for pj > 1 GeV/c; in pp and p-Pb collisions CD- Recently, it has been demonstrated that 
anisotropic flow (V 2 and V 3 ) can be obtained using gluon distributions from classical Yang-Mills simulations of p-Pb 
collisions Ifl2l . Similarly, the AMPT model reproduces the ridge structure in small systems Ifl3l and Pythia 8 Ifffl 
with the tune 4C 03 is able to generate flow-like patterns using color reconnection (CR) and multi-parton interactions 
(MPI) fl 6 ft . It is worth to mention that early LHC results for pp collisions, e.g., the transverse sphericity m and 
the di-hadron correlations ED as a function of multiplicity, revealed the importance of multi-parton interactions to 
describe the data. 

As discussed here na. it is important to look at new observables which allow one to establish the origin of the 
QGP-like effects in small systems. In this context, the present work shows a study of the average transverse momentum 
(p-y) as a function of the hadron mass ( m ) and the event multiplicity (N c h) obtained from different measurements 
reported by experiments at the LHC. Similarities among pp, p-Pb (all multiplicity classes) and peripheral Pb-Pb (60- 
90%) collisions are found, those effects are qualitatively well captured by Pythia 8 , but not by hydro calculations. On 
the other side, for central Pb-Pb collisions the data exhibit a different feature. 

2. Multiplicity dependence of the average pj vs. hadron mass in Pythia 

In Pythia 8 , color reconnection and multi-parton interactions produce flow-like patterns in pp collisions at the 
LHC energies via the creation of boosted strings m. To illustrate this, the MC p-y spectra of different particle species 
(n + + n~, K + + Kr, K^, p + p, 0 -meson, A + A, EL + £ + and LT + D + ) for pp collisions at frs = 7 TeV were studied 
using Pythia 8.201 [T9) (tune Monash 2013 f20j ). The transverse momentum spectra and event multiplicity were 
determined at |y[ < 1 and |/;| < 1, respectively. It has been shown that the Boltzmann-Gibbs blast-wave function HI 
describes within 10 % the p-y distributions of all hadrons under consideration when color reconnection is activated m. 
Actually, the ALICE Collaboration has reported the same level of agreement (within »10%) between the data and the 
blast-wave model for p-Pb collisions GED. Furthermore, with CR the correlation between the parameters obtained 
from the fit to MC p-y spectra, the average transverse velocity, (fry), and the temperature at the kinetic freeze-out, 7f m ; 
behaves like in data HI |23]. Namely, (fry) increases with increasing the multiplicity density (dA' c h/d//), and at the 
same time 7\ m decrease^] Hence, the results from the blast-wave analysis are not enough to draw a final conclusion 
on the possible formation of a strongly interacting QCD medium in high multiplicity pp and p-Pb collisions. 


1 The odd harmonics were not generated because rescattering contributions to the intrinsic correlations were not included. 

2 In simulations without color reconnection, the multiplicity dependence is very weak, i.e., (Jj-y) and 7\ ln do not evolve with multiplicity. 
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Figure 2. Average transverse momentum as a function of the scaled hadron mass ( m/riq ) for minimim bias pp collisions at \fs = 0.2 TeV (left), 
0.9 TeV (middle) and 7 TeV (right) measured by PHENIX and ALICE experiments. 


The main focus of this work is the study of the mass dependence of the average p-\ . In the string hadronization 
model, the string boost is transferred to the final particles when the string breaks up and for the same boost velocity 
a heavy particle will gain more p\ than a lighter one. CR is therefore expected to give a rise of the average pj vs. 
hadron mass. Qualitatively, the effect is in good agreement with data, but so far, the current models of CR are not able 
to reproduce the size of the effect l24j . However, other effects look promising, for example the color ropes formation, 
which by increasing the effective string tension, can also produce flow-like behavior l25l . 

The mass dependence of the average pj has been investigated for two cases, with and without color reconnection. 
When color reconnection is deactivated, the partonic systems fragment independently, and even in this case, the av¬ 
erage pj increases with the hadron mass. It is worth to remember that in LHC data, even for the lowest multiplicity 
event class, where a medium is not expected to be formed, the mass ordering is also seen l26l - l28l . The (pj) increase 
with the hadron mass, therefore is not a distinct feature of hydro, instead it seems to be a characteristic of the fragmen¬ 
tation. Actually, another observation which is reported for the first time is that Pythia 8.201 without CR produces an 
average pj which scales with the reduced hadron mass, i.e., mass divided by the number of quark constituents ( m/n q ), 
and not with the hadron mass: 


m 

(p T ) cc —. (1 ) 

n q 

This is illustrated in Fig. [T] where the correlation between (pj) and m/n q is shown for two event classes: low 
(0 < dA c h/d^ < ( dN ch /dr ])) and high (7 x (dA^h/d^) < dN c h/dp < 8 x (dA c h/d^)) multiplicity. The scaling is valid 
across the multiplicity event classes, though, it is somewhat violated at high multiplicity. Figure [T] also illustrates the 
results when color reconnection is activated, in this case, at low multiplicity the scaling holds, while at high ,V ch it is 
broken and the (pj) for baryons and mesons increase with different slopes. In addition, the mean pj for (f> meson is 
higher to that for protons, hence, neither a scaling with hadron mass would be expected. 

In the hadronization model used in Pythia, the constituent quark number dependence can be built up as follows. 
Quark - (anti)quark pairs are produced during the string break-up process, a transverse momentum is assigned to the 
pair following a flavor-independent Gaussian spectrum. Since no string transverse excitations are assumed, the pj 
is locally compensated between the quark and (anti)quark of the pair. Each hadron is subsequently formed with the 
string segments and the total partonic pj is assigned to the hadron. In this model, a meson corresponds to a small 
piece of string with a quark in one end and an (anti)quark in the other l29ll . This fragmentation is universal, in the 
sense that is the same for all the primary hadrons produced in the collisions, and therefore a dependence of the average 
hadron pj with the constituent quark number is expected to occur. Following the same prescription, in events with 
large number of MPI and with CR, partons from different origin (MPI, initial state radiation, etc.) are allowed to 
recombine to form hadrons. This sort of recombination leads to a general n q dependence of the pj which becomes 
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more complex since, as discussed in the previous paragraph, on top of the universal Gaussian fragmentation the string 
segments are also boosted. In contrast, hydro calculations for small systems Boil give a (pj) which is independent of 
n q . 

The next sections are dedicated to study the same effects in different colliding systems using published LHC data. 

3. Study of (p-\ ) using LHC data 

The scaling of (/? T ) with m/n q is first investigated in inelastic pp collisions at different center of mass energies, 
namely, 0.2 TeV, 0.9 TeV and 7 TeV. Data were obtained from different publications of experiments at the RHIC RD 
and at the LHC f32l]34l . Figure [2] shows the correlation plots, (pj) vs. m/n q , which are obtained. For 0.2 TeV data, 
where the event multiplicity is the smallest, ( dN c h/drj ) ~ 2.25 1351 , the m/n q scaling scaling holds even if heavy 
hadrons like E and Q are considered. The same observation is valid for pp collisions at 0.9 TeV, where the event 
multiplicity is still small, (dN^/dp) ~ 3.81 1361 . For the highest center of mass energy, where the identified hadron 
pj spectra are available, the scaling is slightly broken as observed in Pythia 8.201 with color reconnection. As in the 
model, the (pj) for baryons rises faster with the reduced mass than for mesons. Also the average pj for baryons is 
higher to that for mesons. 

Using the existing ALICE data for p-Pb collisions l26l . the average p-\ for n + + n , K + + K , p +p, K° and 
A + A as a function of m/n q are obtained and plotted in the upper part of Fig. [5] Two multiplicity event classes are 
shown, (dV C h/d/ 7 )= 4.4 and (dV c h/d//}=45. As in Pythia 8 , for low multiplicity events the scaling of (pj) with m/n q is 
observed, while for high multiplicity events it is broken. There, the mean pj for baryons is higher than for mesons. 
The analogous analysis for Pb-Pb data 13814401 is presented in the bottom panel of Fig. [3j there, also the <p meson 
(pj) is reported. Like in Pythia 8 , the scaling of (p-\ ) with m/n q is only valid for the most peripheral Pb-Pb collisions 
(80-90%). For the most central Pb-Pb collisions (0-5%), the mean pj for baryons and mesons as a function of m/n q 
exhibit different trends. In addition, for this event class the average p T for protons is roughly the same to that for <p 
mesons. This observation triggers the question of whether a scaling of (pj) with mass, and not with m/n q , is seen in 
the most central Pb-Pb events. 

The possible universal scaling with m is investigated. Figs. [4] and [5] show the multiplicity dependence of the 
average pj as a function of the hadron mass for p-Pb and Pb-Pb collisions, respectively. The universal mass scaling 
is only observed in the 0-60% centrality classes for Pb-Pb collisions. While, it is broken for the rest of the Pb-Pb 
centrality classes (60-90%) and for all the multiplicity classes measured in p-Pb collisions. In Fig. [4] also hydro 
calculations for p-Pb l30l are shown, there, the mass scaling holds for all the multiplicity classes. 

One thing which can be seen is that for the same multiplicity density, fcLV ch /d/ 7 ) ~35, the slope of (pj) vs m is 
higher in p-Pb than in Pb-Pb collisions. This can be quantified if a first degree polynomial is fitted to data, in this 
exercise the correlations (p T ) vs m/n q for p-Pb and Pb-Pb data were used. Figure[ 6 ]show the results for mesons, where 
the slopes obtained from the fits are plotted as a function of the event multiplicity. For comparison, also the results for 
inelastic pp collisions at different center of mass energies are shown. At low multiplicity, where the universal scaling 
of (pj) vith m/n q is approximately valid, the data points from all the systems follow the same trend, and then, for 
higher multiplicities, p-Pb and Pb-Pb data deviate. 

4. Conclusions 

Collectivity is considered one of the key ingredients to establish the formation of an equilibrated QGP So far, the 
way to verify whether a system behaves collectively or not relies in the success of the hydro models to describe flow 
observables. Namely, the pj distributions of identified hadrons and v n , which are connected with radial and anisotropic 
flow, respectively. In the present work the details of existing pj spectra for different colliding systems were studied 
aiming to check two things. On one side, how precise hydro models describe data, specifically, to check in data the 
mass scaling of the average pj as predicted by hydro Go). And on the other side, to test an alternative explanation 
based on color reconnection and multi-parton interactions. This investigation brought the following results: 

• Low multiplicity pp collisions simulated with Pythia 8 exhibit an universal scaling of (pj) with m/n q (mass 
divided by the number of quark constituents). When CR is switched-off, the effect does not depend on mul¬ 
tiplicity. In the opposite case, this universal scaling is broken, i.e., a deviation from the linear behavior is 
observed, and the average pj for baryons and mesons increase following different lines. 
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Figure 3. (Color online). Average transverse momentum as a function of the hadron mass normalized to the number of constituent quarks. Results 
for the lowest and highest multiplicity event classes in p-Pb collisions are shown in the upper panel. For Pb-Pb collisions, the results for the most 
peripheral and the most central events are shown in the bottom panel. The solid (dashed) line indicates a linear fit to meson (baryon) data. The 
average multiplicity densities, (diV c h/d77), for the different event classes have been measured by ALICE in \tj\ < 0.5 (26l|33. 
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Figure 4. (Color online). Multiplicity dependence of the average transverse momentum as a function of the hadron mass in p-Pb collisions 
at y.s'NN = 5.02 TeV (solid markers). A first degree polynomial (solid line) is fitted to the meson data. For three multiplicity classes, hydro 
calculations (30) are also shown (empty markers), linear fits to meson (pj) calculations are shown as dashed lines. The average multiplicity 
densities, (dA c h/d 77 >, for the different event classes have been measured by ALICE in I 77 I < 0.5 126). 
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Figure 5. (Color online). Centrality dependence of the average transverse momentum as a function of the hadron mass in Pb-Pb collisions at 
-\/^nn = 2.76 TeV. Solid line indicates the first degree polynomial which fits the meson data. The average multiplicity densities, (dA/ch/d^), for the 
different event classes have been measured by ALICE in I 77 I < 0.5 (57). 
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Figure 6. (Color online). For each event class, a first degree polynomial is fitted to the average pj vs. the reduced hadron mass (m/n q ). In this 
figure, the slopes extracted from the fits are plotted as a function of the event multiplicity for different colliding systems. 


• LHC data for p-Pb collisions as a function of the event multiplicity give results which are qualitatively similar to 
those observed in Pythia 8. The same observation holds for inelastic pp collisions at RHIC and LHC energies. 

• For LHC Pb-Pb data the average p-y as a function of mass shows an interesting phenomenon. Up to centralities 
of about 50-60% the baryon and meson (p-y) do not scale with the same dependence (similar to the results 
obtained with Pythia). For more central collisions the scaling with mass seems to work. 

The results imply that the application of the hydro model should be done with great care because the scaling with 
mass is not an universal behavior of the data. The hydro model assumes that the behavior of all particles in the mostly 
pion cloud created after the collision is identical. However, without taking into account the difference in the pion 
- particle interaction cross section (as was observed at SPS and RHIC energies for the slope). Again the results of 
Pythia simulations surprisingly well reproduce the different behavior of mesons and baryons raising the possibility of 
an altogether alternative interpretation. 
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